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Abstract. 

We analyze the suppression of particle production at large transverse momenta in 
(0 — 5% most) central collisions of gold nuclei at */snn = ^00 GeV and lead nuclei 
at ^sp^N = 2.76 TeV. Full next-to-leading order radiative corrections at 0(atg), and 
nuclear effects like shadowing and parton energy loss are included. The parton energy 
loss is implemented in a simple multiple scattering model, where the partons lose an 
energy e = A x dE/dx per collision, where A is their mean free path. We take e = kE 
for a treatment which is suggestive of the Bethc Heitlcr (BH) mechanism of incoherent 
scatterings, e = V aE for LPM mechanism, and e = constant for a mechanism which 
suggests that the rate of energy loss (dE/dx) of the partons is proportional to total 
path length (L) of the parton in the plasma, as the formation time of the radiated 
gluon becomes much larger than L. We find that while the BH mechanism describes 
the nuclear modification factor i?AA ^ or P T — ^ GeV/c (especially at RHIC energy), 
the LPM and more so the constant dE/dx mechanism provides a good description at 
larger px- This confirms the earlier expectation that the energy loss mechanism for 
partons changes from BH to LPM for pr > A < k T >, where A w 1 f m and < k T >sa 1 
GeV 2 is the average transverse kick-squared received by the parton per collision. The 
energy loss per collision at the Js^$ —2.76 TeV is found to be about twice of that at 
0.2 TeV. 



1. Introduction 

The suppressed production of particles having large transverse momenta, known as 
jet-quenching [Tj [2] along with the strong collective flow [3J H] and the success of 
recombination model [5] in explaining the constituent quark number scaling observed for 
the elliptic flow, at the experiments performed at the Relativistic Heavy-Ion Collider 
(RHIC) at Brookhaven National Laboratory have provided a strong evidence for the 
production of quark gluon plasma (QGP) in collisions of nuclei at relativistic energies. 
The jet-quenching [6] and elliptic flow [7] have already been confirmed in collision of 
lead nuclei at A /sj^ r = 2.76 TeV at the Large Hadron Collider. The simplest measure 
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Figure 1. A comparison of inclusive production of neutral pions [IT] with NLO pQCD 
calculations at ^/sjxjjvf = 200 GeV. 



of the jet-quenching at large pt is given by the nuclear modification factor R\\, which 
provides the ratio of particle yield in AA collisions to that in pp, defined as: 

(1/N*£) d 2 N AA /d 2 p T dy 



Raa 



(1) 



iVcoii (l/iV e p v p t ) d 2 N pp /d 2 p T dy 
where N ev ^ is the number of events and iV co n is the number of binary (nucleon-nucleon) 
collisions. In absence of any nuclear modification, -R^A wou ld be unity at large p?. The 
larger density of the medium likely at the higher LHC energies is expected to give rise 
to a larger loss of energy for a given px, which should lead to a larger suppression. Yet 
the expectation that the spectrum of the partons at large pt will fall less steeply at the 
LHC would mean a smaller suppression for the same energy loss of parton of a given 
Pt- Thus indeed the -R^A a ^ LHC shows an interesting behaviour: it decreases with 
increase in in the window of pt ~ 2-5 GeV/ c and then rises to reach a value of about 
0.4 at pt ~ 20 GeV/c. The results for the most central collisions at RHIC energies can 
also be interpreted to suggest the emergence of a similar trend. In the present work we 
suggest that this result is indicative of change of the mechanism for energy loss as pr 
rises beyond 5-7 GeV/c as the Landau Pomeranchuk Midgal effect sets in to suppress 
radiation of gluons. 

In the next section we give the formulation of our treatment and introduce various 
mechanisms of energy-loss used in our work. In section III, we discuss our results at 
RHIC and LHC energies for the 0-5% most central collisions. Finally we give a short 
discussion and conclusions. 
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Figure 2. A comparison of inclusive production of charged hadrons estimated by the 
ALICE collaboration |6] with NLO pQCD calculations at v /sj^" = 2.76 TeV. 



2. Formulation 



2.1. Particle production in pp collisions 

As a first step we check if the next-to-leading order perturbative QCD [8] describes the 
pp data used as a reference in the experimental studies. This is important since, for the 
LHC measurement, the pp data are estimated by using a variety of scaling properties of 
the particle production in hadronic collisions. 

In perturbative QCD, the inclusive cross-section for the production of particle h in 
a nucleon-nucleon collision is given by: 

E h ^-(y/s,p h ) = J dx a J dx b J dzY^Fi(xa,Q 2 )Fj(x b ,Q 2 ) 

D c/h (z,Q})E c ^^ , (2) 

where Fi(x,Q 2 ) is the distribution function for the i-th parton in a nucleon, D c /h is 
hadron fragmentation function at z — Ph/p c , and o~ij^ c x is parton-parton cross-section. 
We include leading order, 0(ot 2 s ) processes, like: 

q + q -> q + q, 

q + q -> q + q, 

q + g -+q + g, 

9 + 9^9 + 9, 
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(3) 



and the next-to-leading order, O(a^), sub- processes such as: 
q + q ->■ q + q + g , 
q + q -> g + g + flf, 

g + ?' -»■ g + </ + 

q + q ^ q' + q' + g, 
9+9^9+9+9 



(4) 



The running coupling constant a s (fi 2 ) 
. 2s 12tt 



calculated at next-to-leading order, is given by 
6(153 -19iV7) In ln(/i 2 /A 2 " 



1 - 



(33 - 2^ ln(// 2 /A 2 ) 



(5) 



(33 - 2N f ) ln(/i 2 /A 2 

where /x is the renormalization scale, Nf is number of flavours, and A = Aqqq. 
We use CTEQ4M structure functions [9] and Binnewies, Kniehl, and Kramer [TO] 
fragmentation functions. We use the factorization, renormalization, and fragmentation 
scales as Q = Pt, though we have checked the results for the particle production in pp 
collisions using scales Q = pr/2 and Q = 2pr as well. 

Our results along with the 7r° data measured by the PHENIX experiment [TT] are 
shown in Fig. [TJ We see a very good description of the experimental results with-out 
any free parameters. Similar results have been repeatedly reported earlier (see e.g. 
Ref. [IT]). 

For the pp centre of mass energy of 2.76 TeV, no experimental data is available and 
the ALICE collaboration has utilized a multi-pronged approach to estimate [6] this using 
several well established scaling relations. This estimate and our calculations are shown 
in Fig. |2j The good agreement of the estimated cross-sections with the NLO pQCD 
results, holds out the hope of getting an accurate measure of the medium modification 
of -R^A in such studies. 



2.2. Medium modification of Rj[j[ 

By now it is very well established that the medium modification of -R^A ^ s ma i nr y 
caused by the energy-loss suffered by partons while they traverse the plasma due to 
collisions and radiation of gluons, before they fragment and the nuclear shadowing. The 
nuclear shadowing accounts for the modification of distribution of quarks and gluons 
inside nuclei compared to that in free nucleons. We use the nuclear shadowing function 
obtained by Eskola, Kolhinen, and Salgado [12] which is known to describe the NMC 
data accurately. 

A substantial body of literature exists on the dynamics of propagation of fast 
partons through the hot and dense medium produced in relativistic heavy ion collisions, 
where they lose part of their energy through collisions and radiation of gluons before 
fragmenting into hadrons. 
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Figure 3. The nuclear modification function for the production of neutral pions for 
the 0-5% most central collision of gold nuclei at the top RHIC energy, e is energy loss 
per collision and the mean free path is taken as 1 fm. The experimental results are 
from the PHENIX collaboration [27]. 
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These developments can be broadly categorized as follows: 

• The Higher Twist Approximation developed by Guo and Wang [13], which is based 
upon the huger twist correction for final state partons in e + A collisions. 

• The so-called AMY formalism developed by Arnold, Moore, and Yaffe [Hj based 
on hard thermal loop resummation in a pertubative plasma. 

• The so-called ASW formalism developed by Armesto, Salgado, and Wiedemann [15] 
which resums multiple soft gluon emission in the BDMPS approach [16] (see later) 
in a finite length medium using Poisson statistics, and 

• The GLV formalism developed by Gyulassy, Levai, and Vitev [17J which considers 
scattering centres in the opacity expansion. 

Several developments tracing the evolution of the entire jet in the medium using 
Monte Carlo jet-quenching modules, e.g. PYQUEN [18], Q-PYTHIA [19J, JEWEL [20], 
YaJEM [21J, and MARTINI [22] have also been reported. Many of these developments 
have been discussed in detailed reviews [231 121] • 

We use a simple phenomenological model [25] to study the evolution of the 
mechanism of energy loss with increasing and centre of mass energy. This has been 
shown to be quite successful in explaining the results for -R^A PS] at RHIC energies. It 
has also been used extensively to illustrate various consequences of jet- quenching, e.g. 
photon-tagged jets and di-hadron correlations etc. 

The advantage of this treatment, as we shall see, lies in giving reasonable description 
of the medium modification in terms of energy loss per collision and the mean-free of 
the parton, which covers an average path length in the medium. We shall see that this 
works quite well and also helps us to distinguish the energy dependence of the energy 
loss mechanism. 

We focus our attention on the central rapidity and central collisions, so that the 
parton propagates in the transverse direction and the azimuthal dependence of the path- 
length etc. can be safely ignored. Taking the path length as L and the mean free-path 
of the parton a as A a , we estimate the probability for a parton to scatter n times before 
escaping as 



Following Ref. [25], we modify the hadronic fragmentation function D c /h(z,Q 2 ) to 
include multiple scattering and the energy loss of the parton in the nuclear medium. 
Assuming that the average energy loss per collision suffered by the parton a is e a , the 
nuclear fragmentation function can be written as: 




(6) 



1 



N 



zD c/h (z,L,Q 2 ) 



Q 



N n=0 



-J2P a (n,L)[z a n D° c/h (z a n ,Q 2 ) 
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where z% = zE^/E^, E® +1 = E® — e°, z^ = p T /e^, n * s the maximum number of 
collisions for which z® < 1, D®, h is the hadronic fragmentation function which gives 
the probability for the fragmentation of a quark or a gluon to fragment into a tt° or a 
charged particle (or any particular hadron in general), and 

N 

C% = J2 P a {n, L) . (8) 

n=0 

Of-course Eft = E^ in an obvious notation. In the Eq. the first term corresponds to 
the fragmentation of the leading parton and the second term comes from the emitted 
gluons having energy e™. The average number of scatterings over the distance L is 
L/X a , where we take L = 1.2(0.5 x ^Vpart) 1 ^ 3 and A a =1 fm. (This value of L should 
suffice for the most central collisions considered here. For large impact-parameters L 
will vary considerably with the azimuthal angle and its average value will be smaller. 
Consequences of this variation along with those of changing e a and A a with quarks/gluons 
will be reported shortly.) As mentioned earlier, we use BKK fragmentation functions 
from ref. [10] and take the factorization, renormalization, and fragmentation scales to 
be equal to pr- 



2.3. Parton energy loss 

We employ three different prescriptions for the the energy loss per collision (e a = 
X a x dE a /dx). The choices are inspired by the excellent discussion in Baier et. al [16J. 
These authors argue that for light quarks and gluons the collisional energy loss (dE/dx) 
is quite small ( 0.2 - 0.3 GeV/fm for quarks and | times this value for gluons atT« 250 
MeV). Thus the energy loss for them is dominated by the mechanism of radiation of 
gluons. 

The radiation of gluons is then best discussed in terms of the formation time of the 
radiated gluon: 

oj , . 

^form ~ T2 ' 

where w > ^ is the energy of the gluon and kx is the transverse momentum of the 
gluon. For small oj (£f orm < A) incoherent radiation takes place over L/X scattering 
centres. This is the so-called Bethe-Heitler regime and one can derive: 

- f „ ^n c \e (10) 

dx it X 

where N c = 3. We shall write e ~ kE for this case and determine k using the -R^A 
measured in the experiment. 

When the formation time if orm is greater than A but less than L, iV con scattering 
centres act as a single source of radiation, where 

^coh = W A (U) 

and 
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Figure 4. The nuclear modification function for the production of charged particles 
for the 0-5% most central collision of lead nuclei at y / sj^"= 2.76 TeV measured at the 
Large Hadron Collider, e is energy loss per collision and the mean free path is taken 
as 1 fm. The experimental data are from the ALICE experiment [5]. 
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Figure 5. Evolution of the mechanism of energy loss from Bethe Heitler to Landau 
Pomeranchuk Midgal as a function of pr at RHIC and LHC. 



where < k\ > co j 1 is the accumulated transverse momentum and ~ ^coh x < &t >■ 
Thus one can write: 



N, 



UJ 



coh 



A < k% > 



(13) 

v t ^ V ^LPM 

where -^LPM = A < fc^ > is the energy which separates the BH regime from the LPM 
regime. Now, one can write: 



dE a s N c 

— ~ ~" -^LPM-^ • 



dx " 7t A V ~L,rm- ■ ( 14 ) 

In the present work, we denote this scenario with e = y/aE and determine a from 
a description of R\\- 
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Figure 6. Doubling of energy loss per collision as y's^N is increased from 200 GeV 
to 2.76 TeV. 



Finally, when the formation time is much larger than the path-length L the 
coherence is complete and the medium induced total energy loss becomes proportional 
to the square of L, and one gets: 



dE _ 

dx 



< k% > 



x -- (15) 

7T A 

We shall denote this case by having e = constant, for a given L. One may need to 
check it for cases having varying L (e.g., centrality dependence of R\pJ- Results for 
such a study would be published shortly. The NLO code of Aurenche et al had earlier 
been modified by Jeon et al [26J to incorporate the medium modified fragmentation 
functions. We use it for the calculations reported here. 



3. Results 

We have already seen that the NLO pQCD calculations accurately describe the results 
for particle production in pp collisions at both the energies under consideration (See 
Figs. US}. 

The nuclear modification factor -R^A obtained in the present work for central 
collision of gold nuclei at t/^nn equal to 200 GeV, using the three energy loss 
mechanisms discussed above are shown in Fig. [3j The so-called Bethe Heitler mechanism 
is seen to accurately describe the nuclear modification for 2.5 < px < 6 GeV/c (top 
panel) . 

It is interesting to note that the condition necessary for the multiple scatterings to 
be incoherent, when the mean free path is about 1 fm requires that the formation time 
of the gluon is less than 1 fm/c. This will be satisfied if the energy of the radiated gluon 
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is less than 5 GeV and < k T >^ 1 (GeV/c) 2 . Thus indeed the partons having energy up 
to 5-7 GeV would populate this px window. The value of k ~ 8%, obtained here is close 
to (but higher than) that obtained by Jeon et al. [26] for 0-10% centrality. However it 
is smaller than the value suggested by the Eq. [TUJ We shall discuss this later. 

The nuclear modification shows a clear change in slope around pt equal to 5 GeV/c. 
We also note that -^lpm ~ 5 GeV as discussed above. Thus it is interesting that the 
energy loss implied by Eq. [H] provides an accurate description of R\\ for 5 GeV/c 
< pr < 12 GeV/c, when we take a as 0.07 GeV 1 / 2 (Fig. [3j middle panel). We again see 
that while the form of the energy loss is given by Eq. [141 the coefficient for the energy 
loss we get is about half as small. 

The statistics for the modification factor is quite poor at larger p?, however the 
average behaviour seen (see Fig. 3, lower panel) is in fair agreement with energy loss 
mechanism given by Eq. [151 though once again the energy loss per collision is smaller 
than that implied by this equation. 

The results for A /sj^ r = 2.76 TeV, (see Fig. @| are also quite interesting. We now 
see (top panel) that the BH mechanism provides a description of the data over a limited 
Pt window of 5 - 7 GeV/c. It is quite likely that at these higher energies the validity 
of hydrodynamics may extend to higher pt than that witnessed at RHIC energies and 
thus the radial flow may affect the pt window up to 5 GeV/c as well as our results. 

The so-named LPM regime is seen to hold over the pt window of 7-14 GeV/c (see 
FigJH middle panel) with a ~ 0.25 GeV 1//2 , which is still smaller (by about 30%) than 
that implied by Eq. [T4l 

The mechanism which gives a constant energy loss per collision is seen to accurately 
describe the data over the largest pt window of 10 - 20 GeV/c (see FigJH lower 
panel). This is interesting as this mechanism is supposed to work for L < L CT where 
L C t = Ay 1 E I '-^lpNi w hich should be true for E > about 125 GeV for L = 5 fm, while 
it seems to be operating even at lower energies. The energy loss per collision is still 
smaller than that implied by Eq. [151 though it is about twice as large as at the top 
RHIC energy. 

In order to bring out the changing mechanism for energy loss with increasing pt 
we now plot the final results for RHIC and LHC energies using the BH and LPM 
mechanisms (see Fig. [51 upper and lower panels). The results at RHIC energies show a 
dramatic change in the mechanism around pt of 5 GeV/c, as indeed one expects from 
the consideration of formation time of the radiated gluons and the mean free path of 
the partons. 

The comparison for the constant energy loss per collision for the RHIC and the LHC 
energies is shown in Fig. [6j We see that while the results at RHIC are only indicative 
of possible emergence of this mechanism at larger pt, the results at LHC clearly favour 
this mechanism for pt > 10 GeV/c. 
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4. Summary and Discussions 

We have adopted a simple and transparent model of energy loss of partons traversing 
the medium created in relativistic collision of heavy nuclei. The fragmentation functions 
of the partons are then modified to account for the energy loss (and also fragmentation 
of radiated gluons) and nuclear modification function -R^A calculated. The energy loss 
per collision is taken to be either proportional to the energy of the parton, square-root 
of the energy of the parton, or a constant, inspired by the treatment of energy loss 
depending on the formation time of the radiated gluon and the mean free path of the 
parton. The reference particle spectra are calculated using NLO pQCD. 

A dramatic change in the energy loss mechanism is seen at pt ~ 5 GeV/c. We find 
that the Bethe- Heitler mechanism of incoherent scattering prevails for p? < 5 GeV/c at 
both RHIC and LHC energies, while the so-called LPM mechanism, which gives energy 
loss per collision as proportional to the square-root of the energy prevails over the pt 
window of 5-15 GeV/c. In an interesting observation we find that the constant energy 
loss per collision suggestive of a total energy loss increasing with the square of the path 
length may already be becoming relevant for pt > 10 GeV/c. This mechanism is seen 
to prevail over the px window of 10-20 GeV/c at y/^NN = 2.76 TeV. 

While the results generally follow the qualitative arguments given in Ref. [16] for the 
form of the energy-loss per collision as the energy of the parton changes, the energy-loss 
coefficients are smaller than those implied in the above work. One obvious difference is 
while the medium produced in these collisions is expanding and cooling the treatment of 
Ref. [16] is basically for a static medium (with static centres of scattering, a la Glauber 
model), though modifications for expansion are available. Obviously, in an expanding 
medium the parton will initially move in a hotter layer and later in cooler layers. The 
energy loss per collision emerging in the present work will then give an average over the 
history of evolution. The expansion of the plasma will also lead to a L larger than the 
radius of the system, at least at LHC energies. However the expanding volume beyond 
the volume of the plasma produced originally will also be cooler and thus the energy loss 
there may be small. It is also likely that the partons produced near the centre would 
have already lost most of their energy by the time they enter this volume, and thus they 
would fragment to low energy hadrons. The boundaries of the different mechanisms 
may also not be quite sharp. 

It will be interesting to see how these mechanisms evolve with the centrality of the 
collision [28] and how do they fare when the dynamics of evolution is accounted for [29|. 
Is this simple but transparent treatment able to account for the azimuthal anisotropy 
of the momentum distribution for non-central collisions at large Pr? These will be 
addressed in a forthcoming publication [30]. We add that, even though the discussion 
of Ref. [16] focuses on radiative energy loss of the partons the treatment in terms of 
energy loss per collision could easily accommodate the collisional energy loss, as long as 
it is small. 

We conclude that the mechanism of energy loss of the partons weaves a rich tapestry 
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depending on the energy of the parton and the properties of the medium. The simple 
model of multiple scattering used in the work is able to reveal the evolution of the 
mechanism of energy loss with px- 
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